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2. Introducción 


Existen varios tipos de interés que ofrecen las entidades financieras, modelizar el compor- 
tamiento de estas es importante para un mejor aprovechamiento de los recursos financieros 
que se requieran, en este material se enfocará en un tipo de tasa de interés como lo es la tasa 
de interés corta instantánea y un modelamiento estocástico de su comportamiento que es el 
modelo de Vasicek y sus variantes. Además de las aplicaciones del modelo en el modelamiento 
de otros elementos como el precio de bonos por ejemplo. 


A lo largo del material se usan algunos términos en inglés para no dar lugar a malas tra- 
ducciones, se omiten las demostraciones de las fórmulas pero se señala la manera de cómo 
calcularlas. Se hará uso de código en OCTAVE-MATLAB sobre todo para generar simulaciones 
de los modelos aplicados y al final se dará solución a problema de aplicación en el que se 
buscan estimar los parámetros del modelo Vasicek dada una muestra de datos. 


2.1. Preliminares 


Las dinámicas de las tasas de interés son equivalentes a los movimientos financieros del 
mercado de bonos. Un bono con cero cupones que paga 1 unidad monetaria al tiempo de 
madurez T', tiene un precio p(t, T') para cada t < T. Las hipótesis básicas sobre este mercado 
son las siguientes : 


= Para cada T' > 0 existe un mercado sin fricciones de T7' bonos. 
“ p(T,D)=1 
= p(t,T) es diferenciable en T 
Luego la tasa corta instantánea al tiempo t está dada por: 
r(t) = f(t,t) 


El proceso de la “cuenta del dinero"se define como: 


B(t) = exp ( A ruda) 


, este proceso puede ser estocástico, y no sólo variable en el tiempo, de hecho es el numerario 
para precios neutrales al riesgo. El precio en el momento O de un valor derivado que paga X 
en el tiempo T es la expectativa de X/B(t), es decir: 


E ep (- ll ruda) x (2.1) 


la expectativa tomada con respecto a la medida de neutralidad al riesgo. En particular, el 
precio en el tiempo 0 de un bono que paga 1 en T' está dado por: 


B(0,T)=E en (- [ was) (2.2) 


Se enfocará principalmente en la dinámica de la tasa corta bajo la medida de neutralidad al 
riesgo. 


3. Modelos básicos y simulación 


El modelo generalizado que describe la dinámica de curva de tasas de interés a través de 
la tasa corta, es el siguiente: 


dr(t) = lr, t)dt + 0(r,t)dW (t) 


Donde: — dr(t):= diferencial o cambio de tasa en un tiempo t 


p(r,t):= tendencia o “drift"de la tasa 
o(r,t):= volatilidad del cambio de la tasa 
dW(t) := movimiento browniano geométrico 
dt := diferencial o cambio del tiempo 
r:= tasa (es función de t) 
t:= tiempo 


3.1. Modelos 


El modelo clásico de Vasicek describe la tasa corta a través de un proceso de Ornstein- 
Uhlenbeck: 
dr(t) = a(b— r(t))dt + odW (t) (Sel) 


, aquí W es un movimiento browniano estándar y tanto a,b y 9 son constantes positivas. 


El drift en 3.1 es positivo si r(t) < b y negativo si r(t) > b; por lo tanto, r(t) es empujado 
hacia el nivel b, una propiedad generalmente referida como reversión a la media. Se puede 
interpretar b como un nivel de tasa de interés a largo plazo y a+, como la velocidad a la que 
se tiende de r(t) hacia b. 


Se considera también el modelo Ho-Lee de tiempo continuo: 
dr(t) = g(tjdt + odW(t) (8.2) 
, con g una función determinista del tiempo. 


Tanto 3.1 como 3.2 definen un proceso gausiano, lo que significa que la distribución conjunta 


de r (t,),..., r (t, ) es normal-multivariante para cualquier t,,...,t,. Ambos definen los procesos 
de Markov y son casos especiales del proceso general de Gauss-Markov especificado por: 
dr(t) = [g(6) + h(8)r(t)] de + o(t)dW (t) (33) 


, con g, h y O todas funciones deterministas del tiempo. 


La ecuación diferencial estocástica 3.3 tiene solución: 


t t 
r(t) = e4Or(0) + 7 AO0HOSg(s)ds + / AH ($) dW (s) 
0 0 
con 


t 
Hi / h(s)ds 
0 
como se puede verificar mediante la aplicación de la fórmula de Itó. 


Debido a que esto produce un proceso gaussiano, la simulación de r (t,),..., 7 (t,) es un caso 
especial de problema general de muestreo de una distribución normal-multivariante. Pero 
es un caso especial lo suficientemente interesante como para merecer consideración. Para 
equilibrar la manejabilidad con la generalidad, el estudio se centrará en el modelo de Vasicek 
3.1 con b variable en el tiempo y en el modelo Ho-Lee 3.2. Ideas similares se aplican al caso 
general 3.3. 


3.2. Simulación 


Para el modelo de Vasicek con b variable en el tiempo, la solución general anterior tiene 
esta forma: 


t t 
r(t) = e *r(0) + a] e —s)ds + o] ed W (s) (3.4) 
0 0 


De manera similar para cualquier 0 < u < t: 


r(t) = 0 (u) +0 il 


u 


t t 
ect s)ds + o] et gy (s) 
De esto se sigue que dado r(u), el valor r(t) se distribuye normalmente con media: 


ei (y) + u(u, t) 
u(u,t) = af, eih(s)ds 


y varianza: 
t B 
o7(u,t) = o] eat gs = E (1 — e2atiu)) (3.6) 
> 20 
Para simular r en los tiempos 0 = ty < t¡ < <t,, se establecen: 
r(tig1) = at (E) + qu (lo, ti41) + 00 (lis ti01) Liga (3.7) 
con Z1,..., Zn  N(0, 1) independientes. 


Este algoritmo es una simulación exacta en el sentido de que la distribución de los r (t,) ,..., r (t,,) 
que produce, es precisamente la del proceso de Vasicek en los tiempos ty, ...,t, para el mismo 
valor de r(0). Por el contrario, el ligeramente más simple esquema de Euler: 


7 (ti+1) = 7 (,) +0: (0 (€:) — y (84) (ligar — ti) + 0 Y tiva — tiLi41 (3.8) 
conlleva algún error de discretización. 


En el caso especial de que b(t) = b, el algoritmo de 3.7 se simplifica a: 


1 
r (si) mo e tig 1=ti) y. (t;) dE (1 q ¿iaa + alo (1 REN e2otigats)) Zi, (3.9) 
Q 


3.3. Versión estacionaria 
Suponga que b(t) =b y a; > 0. Entonces de 3.5: 
E [r(t)] = e “r(0) + (1-e%)b=>b,t >00 
por lo que el proceso r(t) tiene una media límite (valor esperado a largo plazo). También 
tiene una variación límite (vía 3.6) dada por: 
2 2 
[m Var [r(1)|= Jim > (1-4) = 


De hecho, r(t) converge en distribución a una distribución normal con esta media y varianza, 
en el sentido de que para cualquier x € IR: 


P(r(t)<zx)>0 (a) 


con 9 la distribución normal estándar. 


Ko] 00 = o al ns [9] 1) Ra 


hor 
R o 


pp 
w 


3.4. Ejemplo 


Se considera la ecuación 3.9 (modelo de Vasicek) para valores: a+ = 0,1, b= 0,05, 0 = 0,015 
y r(0) =0, con lo cual: 


1 
a (aa) e e0Mti+1=ti) y. (t;) + 0,05 (1 0 Mt) as 00151 (1 NN eA0D(tiit0) Z 4 
con Z1,...,2n  N(0, 1) independientes , a continuación se presenta la gráfica de 15 simula- 
ciones (figura 1) realizadas por medio del código 1. 


Figura 1: Simulaciones del modelo de Vasicek 
Además de las simulaciones se muestra el valor del interés a largo plazo: b=0.05 
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Código 1: Simulaciones del modelo de Vasicek 


clear ; close all; clc 

pkg load statistics 

delta_tiempo = 0.01;/% dt 
tiempo_final = 20; 

tiempo = O:delta_tiempo:tiempo_final; 
alpha = 0.1; Z Q 

sigma = 0.015; % 0 

19 = (0068 4 lo 

OO AE NO 


soja ls E asilo 


T(1) = tiempo(1); 
X(1) = x0; 
LAZO 
for i = 2:length(tiempo) 
T(i) = tiempo(i); 
Z(i) = normrnd(0, 1); 
X(i) = exp(-alpha+* (tiempo (i)-tiempo(i-1)))*X(i-1)... 


+bx* (1-exp(-alpha+ (tiempo(i)-tiempo(i-1))))... 
+sigma*sqrt ((1/(2*alpha))*(1-exp(-2*alphas (tiempo (1) -tiempo(i-1)))))*Z(1); 
end 
ir) 29) 
hold on; 
end 
plot([0,tiempo_final], [0.05,0.05], "-r"); 
xlabel("t");ylabel("r(t)=1£(t,t)");text([0.01], [0.052], "b=0.05");axis("square"); 


4. Precios de bonos 


Como ya se señaló, los parámetros drift dependientes del tiempo se utilizan típicamente 
para hacer un modelo de tasa corta consistente con un conjunto observado de precios de 
bonos. La implementación del algoritmo de simulación 3.7 está, por tanto, vinculada a la 
calibración del modelo a través de la elección de la función b(t). Lo mismo se aplica a la 
función g(t) en el modelo de Ho-Lee y como este caso es un poco más simple, se considera 
primero. 


El punto de partida es la fórmula de fijación de precios de los bonos 2.2. La integral de r(u) 
de 0 a T que aparece en esa fórmula se distribuye normalmente porque r(u) es un proceso 
gaussiano. De ello se deduce que el precio del bono es la expectativa de el exponencial de 
una variable aleatoria normal. Para una variable aleatoria normal X = N (m,u?), se tiene 
E [exp (X)] = exp (m + (v?/2)), entonces: 


Beso (-/ arias) = exp (-8 || ra +Var |/ ria) (41) 


Para encontrar el precio del bono, por lo tanto, se necesita encontrar la media y la varianza 
de la integral de la tasa corta. 


En el modelo de Ho-Lee, la tasa corta viene dada por: 


r(t) = r(0) + / g(s)ds + o W (t) 


y su integral por: 


[rw =r(0)T + 1 JU strasda + AÑ W(u)du 


Esta integral tiene media: 
T u 
nor+/ J storasda 
o Jo 


Var [0 fa W(w)du! =. Ze. Peal hs Cov [W (u), W(t)] dudt 


0 
= 20? e ududt (4.2) 
ZA 
ss e 


y varianza 


Sustituyendo estas expresiones en 4.1, se obtiene: 


E [exp (- fo r(ydu)| 


exp (=r(0)7 _ E Lo g(s)dsdu + ZP) 


E (0, TE) 


Si se da un conjunto de precios de bonos B(0,T') en el tiempo 0, el objetivo es elegir la 
función gy para que se mantenga esta ecuación. Para llevar a cabo esto se escribe: 


B(0,T) =exp (- í (0.0) 


con f(0,t) la tasa forward instantánea para el tiempo t en el tiempo 0. La curva directa 
inicial £(0, T') captura la misma información que los precios iniciales de los bonos. Igualando 
las dos expresiones para B(0,T”) y tomando logaritmos, se tiene que: 


nor+ so E=/ sona 


Diferenciando dos veces con respecto al argumento de madurez 11”: 


0) 
OT 
Por lo tanto, los precios de los bonos producidos por el modelo de Ho-Lee coincidirán con un 


conjunto dado de precios de los bonos B(0,T') si la función g está vinculada a la curva inicial 


a plazo f(0,T'). Así: 


g(t) = F(0,T) lino: +0% (4.3) 


dr(t)= E (0, T) lr: +0 +) dt + 0dW(t) (4.4) 


Se puede realizar una simulación genérica del modelo Ho-Lee con función drift g escrita como: 
ti41 
YT (ti+1) = YT (t;) + J g(s)ds +0 biy1 = tiZiy1 
t; 


con Z1,Z2,... variables aleatorias independientes N(0, 1). Con g elegido como en 4.3, esto se 
simplifica: 


2 
O 
Y (tig1) = 7 (15) + [F(0, tiz1) — $(0,t0)] + [tia E] + 0 tigo — ti Zig 


Por tanto, no es necesaria ninguna integración de la función drift g; para ponerlo de esta 
manera, cualquier integración que sea necesaria ya debe haber sido tratada en elegir la curva 
a plazo f(0,t) para que coincida con un conjunto de precios de bonos. 


La situación es aún más simple si se requiere que la tasa corta simulada sea consistente solo 
con bonos que vencen en los tiempos de simulación t;,...,t,. Para satisfacer este requisito se 
puede debilitar 4.3 a la condición: 


2 


0 l g(s)ds = F(0,ti41) =F (0,t:) + 5 [A ti] 


Excepto por esta restricción, la elección de g es indiferente; se puede llevarla a ser continuo 
y lineal por partes, por ejemplo. De hecho, ni siquiera se necesita especificar y porque solo 
su integral sobre los intervalos (t;,t; + 1) influye en los valores de r en la malla de tiempo 
ti, tn. 


4.1. Bonos en el modelo Vasicek 


Una solución similar, aunque menos explícita, se aplica al modelo de Vasicek. La integral 
de la tasa corta se distribuye normalmente de nuevo; hay que encontrar la media y varianza 
de esta integral para encontrar el precio de un bono usando 4.1. Usando 3.4, para la media: 


E | fa ria] SE [r(t)] at 
Ez Cd (1 = al 0O+a jo fe s)dsdt 


(4.5) 


9 


Para la varianza se tiene: 


Var | Pr pe =2 Ñ ' Cov (r(t), r(1)] dudt (4.5) 


De 3.4 y para u < t: 


Covlritriw] = [eos 
E a (pat t)_ e lu— ») (4.7) 


Integrando dos veces más como se requiere para 4.6 da: 


Var fra] = a 7 + a (1 e) + . (e —1) (4.8) 


Combinando 4.5 y 4.8 como en 4.1, se llega a una expresión para el precio del bono B(0,T). 


Se observa que 4.8 no depende de r(0) y 4.5 es una transformación lineal de r(0). Luego: 
At,T)=2 (1-0) 
; Q 


CMT = ee e e b(s)dsdu 
rs) y 


, esta ecuación puede simplificarse de esta manera: 


2 | 


(esta 1) 


C(t,T) = (A(t,T) - (T —t)) (s - 2) - ca DP 
, luego sustituyendo 4.5 y 4.8 en 4.1 produce: 
B(0,T) =exp(-—A(0, 7) r(0) +C (0, T)) 
De hecho, los mismos cálculos muestran que 
B(t,T) =exp(-A(t DT) r(t) +C(t, T)) (4.9) 
En particular, log B(t, T') es una transformación lineal de r(t). 


Como en la discusión del modelo de Ho-Lee, la función b(s) se puede elegir para igualar un 
conjunto de precios B(0,T') indexados por T. Si solo interesa igualar un conjunto finito de 
precios de bonos B(0, t,),..., B(0,t,,), entonces solo los valores de las integrales 


ti41 
/ ti) (s)ds 
ti 


necesita ser especificado. Estos son precisamente los términos p(t;,t, + 1) necesarios en el 
algoritmo de simulación 3.7. Por tanto, estas integrales son subproductos del ajuste el mo- 
delo a una estructura temporal y no se requieren cálculos adicionales únicamente para la 
simulación. 


10 


4.2. Simulación conjunta con el valor de descuento 


La mayoría de las aplicaciones que requieren la simulación de un proceso de velocidad 
corta r(t) también requieren valores del factor de descuento 


50 be (- [ wa) 


o, equivalentemente de 


Dados los valores r(0),r (t,),...,7 (t,) de la tasa corta, uno puede, por supuesto, generar 
valores aproximados de Y (t;) usando 


1 


y (E5-1) [£; E tii .to=0 


j=i 


o alguna otra aproximación a la integral de tiempo. Pero en un modelo gaussiano, el par 
(r(t), Y (t)) son conjuntamente gaussianos y a menudo es posible simular caminos del par sin 
error de discretización. Para llevar a cabo esto simplemente hay que encontrar las medias, 
las varianzas y la covarianza de los incrementos de r(t) e Y (t). 


Se ha determinado (ver 3.7) que, dado r(t;,) 
r(ti41) e N (ota (E) + pa (ta, ti41) , 0% (ts, €:31)) 
A partir de los mismos cálculos usados en 4.5 y 4.8, se tiene que, dado r(t;) y Y (t;): 
Y (tig1) e N (Y (€) + py (ts, 141) , 0% (ti, t:41)) 


con 


1 ¿+1 
HY (t;, ti41) == (1 = eli ts) )+ ira f f* iS E )dsdu 
Q 


a? 1 a Dr 
o (ti, ti41) = a (ta =t)+ > (1 — e ?alti+rid) y E Ca 1)) 


Solo resta determinar la covarianza condicional entre r(t;+1) e Y (t;41) dado (r(t;), Y (t;)). 
Para ello se procede de la siguiente manera: 


Cov [r(t), Y (t)] = fy Con tr (1), WSL 


= 2 fe eau Ydu 
= 23 [1 + ne E qa. 


La covarianza requerida viene dada por: 


2 


La correlación correspondiente es 


dv esti) 
O litus)Jor llo ba) 


Pry (li, big1) = 


11 


Con esta notación, el par (r, Y) se puede simular en los tiempos t;,..., t, sin error de discre- 
tización utilizando el siguiente algoritmo: 


rías). = Hr (2) + platas) +0 (ita) A (1) 
Y (ti41) Y (8) + py (ti, ti+1) + 0Y (ti, €i41) [Pry (to, £i41) Zi(1 +1) 
+y 1 pry (ti, ti41)Za(i +1) 


donde (%,(3), Z2(1)),1 = 1,...,n, son bivariables independientes normales estándar vectores 
aleatorios. 


4.3. Cambio de numerario 


Hasta ahora, se ha considerado la dinámica de la tasa corta r(t) solo bajo la medida 
de neutralidad al riesgo. El activo numerario asociado con el medida neutral al riesgo es 
P(t) = exp e r(u)du y la característica definitoria de este medida de probabilidad es que 
hace que los precios descontados de los bonos sean B(t, T)/8(t) martingalas. De hecho, para 
la dinámica de los precios de los bonos bajo el régimen gaussiano se consideran modelos de 


la forma (para T fijo): 


dB(t,T) 
o r(t)dt — A(t, T)odW (t) (4.10) 


con A(t, T) determinístico; esto se sigue de 4.9. La solución de esta ecuación es: 


t 1 t 
B(t,T) = B(0,T)exp (/ (0) =- JA Td - o] Alu, TW (4) 
0 0 
de lo cual es evidente que 


B(t,T) 
B(t) 


es una martingala exponencial. 


= B(0, T) exp (-37 / A?*(u, T)du — o [| atu. rjav() (4.11) 


La “forward measure"para cualquier fecha Tp es la medida asociada con la toma del bono 
B(t,Tp) (Tp —vencimiento) como numerario activo. La característica definitoria de la medida 
a plazo es que hace que las razones B(t, T)/B(t, Tp) martingalas para T < Tp. Está definido 
por la razón del proceso de verosimilitud: 


(7) B (t, Tr) B(0) 


dPs ),  B)B(0,Tp) 


y esto se da en 4.11 hasta un factor de 1/B(0,T). Del teorema de Girsanov's, se deduce que 
el proceso W** definido por 


dW*E (t) = dW(t) +04 (t,Tp) dt 


es un movimiento browniano estándar bajo Pr,. En consecuencia, la dinámica de la Modelo 
Vasicek se convierte 


dr(t) = a(o(t) — r(t)) dt + odW(t) 
t)) dt + 0 (AW* (t) — 0 A (t, Tp) dt) (4.12) 


= alb(t) - A (t,Tp) — r(t)) de + odWT"(t) 


12 


Por lo tanto, según la medida a plazo, el proceso de tipos cortos sigue siendo un proceso de 
Vasicek pero el nivel de reversión b(t) se convierte en b(t) — a?2A(t, Tp). 


El proceso en 4.12 se puede simular usando 3.7 con b(t) reemplazado por b(t) — a?A(t, Tp). 
En particular, simulando W** de la misma forma que se hace con cualquier otro movimiento 
browniano estándar. El algoritmo de simulación no “sabe"que está simulando un movimiento 
browniano bajo la medida de avance en lugar de bajo la medida de neutralidad al riesgo. 


Suponiendo que se quiere fijar el precio de un valor derivado obteniendo un pago de g (r (Ty»)) 
en el momento Tp. Bajo la medida de neutralidad al riesgo, se fija el precio del valor compu- 


tando a 
E [elo * 029 (r (Tp))] 


De hecho, g podría ser una función de la ruta de r(t) en lugar de solo su terminal valor. 
Cambiando a la medida de avance, esto se convierte en 


T m 
Er, e” ñ F r(1)du y (r (Tp)) (AL), | Er E So F r(u)du y (r (Tp) ee 
F 


eN 
B (0, Tp) En, [9 (r (Tr) 


4.4. Ejemplo 
Dadas las ecuaciones la ecuación de B(t, T') (ver 4.9), se tiene para a = 0,5,b = 0,06, 0 = 
0,02,T =2 y ro =0,02: 
B(t, 2) = exp (—A (t, 2) r(t) + C (t, 2)) 


donde: 
1 


0,25 


A) == =P 


C(t,2) =(4(t,2) - (2-1)) (0.00 _ 002 ] 02 


20,25) (0,25) 
, a continuación se muestran las gráficas de 15 simulaciones para el precio de los bonos con 
estos parámetros (figura 2), esto por medio del código 2. 


A(t, 2) 
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Figura 2: Simulaciones del modelo de precio de bonos 
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Código 2: Simulaciones del modelo de precio de bonos 


clear close all cle 


pkg load statistics 


delta_tiempo = 0.01;f% dt 
tiempo_final = 2; 


tiempo = O:delta_tiempo:tiempo_final; 


alpha = 0.25; Z Q 
sigma = 0.02; / O 
PISIOSOS 40D 

XO = 0.02; 4 ro 


oe sE alla 
T(1) = tiempo(1); 
X(1) = XO; 
ENERO 


A(1) = (1/alpha)*(1-exp(-alpha*+tiempo_final)); 


C(1) = (A(1)-tiempo_final)*(b-(sigma”2)/(2*alpha”2))... 


-((sigma”2) / (4ralpha))*A(1)72;; 
B(1) = exp A(1)*X(1)+C(1)); 


for i = 2:length(tiempo) 


14 


T(i) = tiempo(i); 

Z(i) = normrnd(0, 1); 

X(i) = exp(-alpha* (tiempo (i)-tiempo(i-1)))*X(i-1)... 

+bx* (1-exp(-alpha+ (tiempo(i)-tiempo(i-1))))... 

+sigma*sqrt ((1/(2*alpha) )*(1-exp(-2*alpha* (tiempo (1) -tiempo(i-1)))))*Z(1); 


A(ií) = (1-exp(-alpha*(tiempo_final-tiempo(i))))/alpha; 
C(i) = (A(i)-(tiempo_final-tiempo(i)))*(b-(sigma”2)/(2*alpha”2))... 
-((sigma”2) / (4ralpha))*A(1)72; 


B(1) = exp CA(D*X()+C(4)); 
end 
plot(T, B); 
hold on; 
end 
xlabel("t"); 
ylabel("B(t,T)=exp(-A(t,TD)r(t)-C(t,T))"); 


axis("square"); 
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5. Problema de aplicación 


5.1. Planteamiento 


Se consideran los siguientes 20 valores que corresponden a intereses registrados durante 5 


meses por periodos homogéneos de tiempo: 


i 1 2 3 4 9 6 T 8 9 10 
Interés | 3 | 1.2693 | 1.196 | 0.9468 0.9532 | 0.6252 | 0.8604 | 1.0984 | 1.431 | 1.3019 

i 11 12 13 14 15 16 17 18 19 20 
Interés | 1.4005 | 1.2686 | 0.7147 | 0.9237 | 0.7297 | 0.7105 0.8683 | 0.7406 | 0.7314 73140 


Se desean generar 50 simulaciones que muestren comportamiento del interés si se suponen 


las mismas condiciones para los próximos 5 meses por medio del modelo Vasicek 


5.2. Solución 


Se considera la ecuación 3.9: 


1 
A (besa) — e7Ati41—ti) yo (t,) +b (1 e O ae az (1 _ e2oltigi=ti)) Z yy 
Q 


Sean: 0t = tiy1 — ti, 141 =P (0; +1), 1, = r(t,) y Z = Z;,1, luego la ecuación se cambia a la 


forma: 


1 
ra = e + b (1 po dla af (1 _ aid Al 
Q 


La relación entre los r,+1 y r, es lineal con una variable aleatoria normal añadida, es decir: 


Ti = Ar +B+e 


2a 


, donde 
A = et . B=b (1 e eaót) E 1—e-2a8t 
, por lo cual: 
A a A) 


Luego usando el ajuste de mínimos cuadrados, se tiene que: 


o (riera) > - (ica Ti Dir 1441) 
( 


A 
Pa 1 me a a o ¿ Ts) 


Se realizan los cálculos y las simulaciones con el código 3: 
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Código 3: Estimación de parámetros para el modelo Vasicek dada una muestra 


clear 


close nal cite 


pkg load statistics 


R 
ll 


n 


length(r)-1; 


[3 1.2693 1.196 0.9468 0.9532 0.6252 0.8604 1.0984 1.431]; 
[r 1.3019 1.4005 1.2686 0.7147 0.9237 0.7297 0.7105 0.8683 


delta_tiempo=0.25; 


X=r (1 
Y=r (2 


YX 


5 


BAS 


RE 
<< 
< 

T 


n); 
:end) ; 


sum(X) ; 
sum(Y) ; 
sum(X.*X) 
sum(X.*Y) 
sum(Y.*Y) 


A = (n*rxy-rx*ry)/(n*rxx-1x72); 
B = (ry-Axrx)/n; 


e = sqrt((n*ryy-ry"2-(A*rxy-rx*ry))/(n*(n-2))); 


alpha 
b=B 


sigma 


tiemp 


= -log(A)/delta_tiempo; 


CAD: 


= exsqrt(-2*1o0g(4) / (delta_tiempo*(1-A72))); 


o_final = 


n*delta_tiempo; 


tiempo = O:delta_tiempo:tiempo_final; 


subpl 


OZ) 


plot (tiempo,r); 


hold 


on; 


legend ("muestra"); 


G=[ 
Ol 


o ls 


= 1:length(G) 
for k = 1:G(1) 


tiempo(1); 


TAN 
X(D =r(1); 
ZNASZRO: 


for i = 2:length(tiempo) 


T(i) 
Z(i) 
xG) 


tiempo(i); 


normrnd(0, 1); 


exp(-alpha+ (tiempo (1) -tiempo (1-1)))*X(i-1)... 
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0.7406 0.7314 0.6232]; 


+bx* (1-exp(-alpha+* (tiempo (i)-tiempo(i-1))))... 
+sigma*sqrt ((1/(2*a1pha) )* (1-exp(-2*a1pha* (tiempo (1) -tiempo(i-1)))))*Z(i); 


end 
subplot(1,2,1) 
plot 
hold on; 
alabelUEINA 
ylabel (Mr 
axis("square"); 

end 

subplot(1,2,1); 


end 


parámetros 
alpha 
b 


sigma 


hesperanza teórica 


EE = ((7(1)-b) .*exp((-alpha) .*+tiempo)).+b; 


fvarianza teórica 


pvar = (sigma”2)/(2*alpha) ; 


var = (-pvar).*+exp((-2*alpha).*+tiempo) .+pvar; 


plot(tiempo, EE,"-k" ,'LineWidth',2); 


plot(tiempo, EE.+(2.*sqrt(var)), "--k" ,'LineWidth',2); 
plot(tiempo, EE.-(2.*sqrt(var)), "--k" ,'LineWidth',2); 


El código 3 dejó los siguientes resultados: 
a=5,11617 b=0,9206  0=3,7245 


con lo cual la ecuación del modelo Vasicek teórico para r(t) queda así: 


(1 e 2(5,1617)(0,25)) 7 
51617 | i ) 


ri = e 76416170025), + 0,9206 (1 A) ep 27245 > 
, además el código 3 también muestra dos imágenes (ver figura 3), en la izquierda están las 
gráficas de la muestra real y la gráfica de una simulación del modelo teórico, en la izquierda 
aparecen las gráficas de 50 simulaciones del modelo y la gráfica de la esperanza E(r(t)) junto 
con las gráficas E(r(t)) + 2yVar(r(t)) 
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Figura 3: Estimación de parámetros del modelo de Vasicek dada una muestra 


3 T T T T 6 T T T T 
—— Muestra 
— data1 


ION de 
IN 0 


5.3. Análisis de resultados 


= Se nota en las 50 simulaciones que los cambios en el valor del interés varía bastante de 
forma repentina 


= Teóricamente se espera que el valor del interés se acerque a la curva de esperanza, esto 
se logrará en un periodo más largo de tiempo 


= Se nota que las curvas de esperanza y esperanza mas y menos varianza son representa- 
tivas del comportamiento del interés en este corto periodo de tiempo 
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6. 


Conclusiones 


En el modelo de Vasicek se observó el fenómeno de reversión a la media, que ayuda a 
predecir los movimientos futuros de las tasas de interés 


Cuando el factor drift (a(b— r(t))) es menor que 0, es probable que la tasa de interés 
disminuya, por el contrario si el factor drift es mayor que 0, es probable que la tasa de 
interés aumente 


Para que el modelo de Vasicek funcione de manera estable, el parámetro a (es decir, la 
velocidad de reversión media) siempre debe ser positivo. 


La volatilidad del mercado (o riesgo de mercado) es el único factor que afecta los 
cambios en las tasas de interés en el modelo de Vasicek (utiliza parámetros fijos). Sin 
embargo, múltiples factores pueden afectar la tasa de interés en el mundo real, lo que 
hace que el modelo sea menos práctico 


Para el desarrollo del modelo Vasicek y sus derivados ha sido esencial el denominado 
Cálculo de Itó, y en particular, de su integral estocástica en términos de la cual se 
expresa la solución del modelo de Vasicek. 


El modelamiento para el precio de bonos supone poder generar también una curva de 
rendimiento de bonos conocida como yield curve 
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